Fermi Velocity Enhancement in Monolayer and Bilayer Graphene 
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In single- layer graphene sheets non-local interband exchange leads to a renormalized Fermi-surface 
effective mass which vanishes in the low carrier-density limit. We report on a comparative study of 
Fermi surface effective mass renormalization in single-layer and AB-stacked bilayer graphene. We 
explain why the mass does not approach zero in the bilayer case, although its value is still strongly 
suppressed. 



I. INTRODUCTION 

Graphene is an atomically thin two-dimensional elec- 
tron system composed of carbon atoms on a honeycomb 
lattice. Experimenters have recently^ made progress with 
techniques which enable isolation and study of systems 
with one or a small number of graphene layers. Much 
of the fundamental physics interest in graphene sys- 
tems follows from the fact that its envelope-function low- 
energy Schrodinger equation is equivalent to the massless 
limit of a two-dimensional Dirac equation. In the case 
of graphene the spinor structure in the Dirac equation 
refers to honeycomb-sublattice and Brillouin-zone val- 
ley, instead of spin and electron-positron, degrees of free- 
dom. Graphene therefore presents a new type of many- 
body problem in which single-particle quantum mechan- 
ics is relativisitc, but interactions are essentially non- 
relativistic instantaneous Coulomb interactions. 

One of the key consequences of electron-electron in- 
teractions in single-layer graphene is an enhancement of 
the quasiparticle Fermi velocity^ii^ which diverges loga- 
rithmically as the carrier density goes to zero. For mod- 
erate interaction strengths^i^ the velocity enhancement 
(Fermi liquid effective mass suppression) is dominated by 
non-local interband exchange effects which are most eas- 
ily explained using the pseudospin language'^'^ explained 
below. In this article we present a close comparison be- 
tween the physics of mass renormalization in single-layer 
and AB-stacked bilayer graphene, explaining why the ve- 
locity renormalization is finite in the low-density limit in 
the latter case. 



II. PSEUDOSPIN DESCRIPTION OF 
VELOCITY RENORMALIZATION IN 
SINGLE-LAYER GRAPHENE 

We begin by recalling the continuum-model Hamilto- 
nian {Ti — 1) of massless Dirac fermions interacting with 
non-rclativistic Coulomb interactions. The Hamiltonian 
has one-particle band-energy and two-particle Coulomb- 
interaction contributions: 
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where v is the bare Fermi velocity, a, (3 are labels for the 
two honeycomb sublattices which we treat as the labels of 
a quantum spin-l/2-like (pseudospin) degree-of-freedom, 
S is the area of the system, and Vq is the 2D Fourier 
transform of the interparticle potential. (Because our 
analysis of Fermi velocity renormalization is based on 
a screened exchange approximation, in which different 
spins and valleys are independent, we do not explicitly 
account for these degrees of freedom.) From Eq. H]) we 
see that in pseudospin language the band Hamiltonian 
consists of a momentum-dependent pseudospin effective 
magnetic field which acts in the direction of momentum 
fc. The band eigenstates in the positive and negative 
energy bands have their pseudospins either aligned or 
opposed to the direction of momentum. 

When interactions are treated in a mean-field approx- 
imation we obtaio^ 



(2) 
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where the density matrix pap{k') = (^'olV'fc ^■^fc.al^'o) 
and l^'o) is the mean-field-theory ground state. (In both 
single-layer and bilayer cases we assume that transla- 
tional symmetry is not broken so that the electron density 
is constant and the Hartree contribution to the mean- 
field Hamiltonian can be ignored.) We parametrize the 
pseudospin-density matrix pai3{k) in terms of charge and 
pseudospin-density contributions : 



Pa/3(fc) = 



Sal: 



AO) 



n{k) ■ crpa 



(3) 



where n^^), are noninteracting band occupation factors. 
Assumming that the valence band (label s = —1) is full 
and that spatial isotropy is not broken, the occupation 
factors are completely fixed by the carrier density n. For 
example, for a n-doped system at T = rS^^ = Q{k-p — k) 



k,a,f3 
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for s 



'1 and n 



(0) 
k.s 



1 for s = — 1, where kp oc is 
the Fermi wavevector and 0(x) is a step function. The 
unit vector h{k) specifies the pseudospin orientation of 
the positive energy band at momentum k. 

Using Eq. ([3]) in Eq. ^ and adding the band ki- 
netic energy one can easily find that the total mean-field 
Hamiltonian is 

Hhf - ^ij^c,f3Bo{k) + (Tap ■ B{k)]i^k.p, (4) 
where the self-consistent Hartree-Fock fields are defined 

by 



Bo{k) 



(Pk' 



vk-k' - 1/2] 



and 



B{k)^vk- j ^vk-k'f-{k')n{k') 



(5) 



(6) 



In Eq. (O we have subtracted the neutral system value 
of Bq, which is easily seen to be wavevector independent. 
This subtraction is necessary in order to make Bq finite 
and represents only a convenient choice for the zero of 
energy. In Eq. ^ the first term is the band-structure 
pseudospin magnetic field, and the second term is the 
exchange field. We also have introduced the short-hand 
notation 



(0) , (0) 



(7) 



which recognizes that n^^\ depends only on k 

We now demonstrate that the solution of the Hartree- 
Fock equation is 



Ifcl 



B{k) =vk- 



<Pk' 
(2^ 



Vk-k'f-{k')k' 



(8) 



This result is a consequence of the isotropic interaction 
Vk-k' = v{\k — k'\) which can be decomposed in angular 
momentum components using 
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(9) 
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with 



X v{^/k^ + k'^ -2kk' cos{9)) . (10) 

Substituting this representation in Eq. ([8|), performing 
the angular integration over ipki we find that 



'jk — 



dk' 
2^ 



k'f^{k')V,{k,k') 



k . (11) 



The ultraviolet cut-off fcmax is required^ because the in- 
tegral is logarithmically divergent at large k' . We set 
fcmax to a value ~ 1/a, where a ^ 0.246 nm is graphene's 
lattice constant, set by the energy scale over which the 
Dirac model is valid. As expected S°'^(fc) is oriented 
along k and \B'"^{k)\ — vk + S(fc) depends only on k. 
Here we have introduced the pseudospin-dependent part 
of the Hartree-Fock self-energy 



S(fc) 



"""^ k'f_{k')V,{k,k') . (12) 
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Note that the fact that quasiparticles in monolayer 
graphene have chirality J = 1 has selected the first mo- 
ment Vi{k,k') of the Coulomb interaction. 

For convenience we assume a n-doped graphene layer 
and take T = so that 

A(fc) = i[e(fcF-A:)-i] - -le(fc-fcF) 

f+{k)- 1/2 - ^eikp-k) (13) 

where kp is the conduction-band Fermi radius. The 
renormalized conduction s = -I- and valence s = — quasi- 
particle energies are therefore given by 

Es{k) ^svk-— [ ^ dk'k'Voik, k') + sE(fc). (14) 
471- Jo 



with 



S(fc) 
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dk'k'Vi{k,k'). 



(15) 



The first interaction correction in Eq. (fT4|) is similar to 
the familiar self-energy correction which occurs in an or- 
dinary two-dimensional electron gas. When the Coulomb 
interaction is not screened it (famously and incorrectly) 
predicts a divergence of the quasiparticle velocity at the 
Fermi energy. Our focus here is on the second self-energy 
correction, exhibited explicitly in Eq. (|15p . which is re- 
sponsible for enhanced splitting between valence and con- 
duction bands. To account partially for screening and 
to avoid the well-known Fermi velocity artifact of mean- 
field-theory in systems with long-range interactions we 
have used a screened interaction potential of the form 



Vk 



e{XqTF + k) 



(16) 



where qxF = Qfccfcr is the Thomas-Fermi screening vec- 
tor and A is a control parameter which can be adjusted 
to represent the weaker screening at higher frequencies, 
(aee = e^/efiv is the graphene continuum- model fine 
structure constant.) We will mainly be interested in the 
contribution to the self-energy at fc ~ fcp due to exchange 
interactions with states at fc' ^ fcp, for which A plays lit- 
tle role. 

The dependence of S(A:) on k/kp for aee = 0.5 and 
A = kynax/kp = 10 is illustrated in Fig.[T] In this plot we 
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FIG. 1; (Color online) The pseudospin-dependent part of the 
Hartree-Fock self-energy E(fc) (in units of ep) for monolayer 
graphene as a function of fc/fcp (for wave vectors up to the 
cut-off) for ftcc = 0.5 and A = 10. The vertical (black) line 
indicates the point k = kp. The data labeled by the dashed 
(red) line correspond to A = 10~^ in Eq. (|16p . while the data 
labeled by the solid (blue) line correspond to A = 1. Note the 
rapid change of E(fc) around fc = fcp when A = 10"^ which is 
an artifact of theories which do not accout for screening. The 
decrease in E for k near fcmax is an artifact of the ultraviolet 
cut-off, but plays little role in the region of interest fc ~ fcp. 



have shown S(fc) for two different values of the control 
parameter A in Eq. p6p : A = 10~^, corresponding to an 
essentially unscreened Coulomb potential, and A = 1 cor- 
responding to the Thomas-Fermi screened potential. We 
see that the unphysical A = anomaly at fc = fcp emerges 
only at very small values of A. For non-zero values of A, 
I](fc) is a smooth function of k and its contribution to 
the quasiparticle velocity at fc = /cp can be obtained by 
expanding S(fc) for k ^ fcmax- For A = (unscreened 
Coulomb interactions) this expansion^ can be performed 
analytically. We find that for x > x' 

V„i{xkF,x'kF) ^ — — Vm{x,x') (17) 
efcp 

where 

/"OO 

Vm(x,x') = / dt J,n{tx)Jm{tx') 

Jo 

x"" r(m-fl/2) 

^.m+i r(m + i)r(i/2) 

X 2Fi{m + l/2,l/2,m + l,x'^ /x"^). (18) 

Here Jm{x) is a Bessel function of order m, V{z) is the 
Gamma function, and 2Fi{a^b,c^ z) is the Gauss hyper- 
geometric function. (The x < x' result is obtained by 
interchange x x' .) For k^ «C fcmaxj a requirement for 
the utility of the Dirac continuum model of graphene. 
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FIG. 2: (Color online) Pseudospin independent exchange en- 
ergy contribution Bo{k) (in units of ep) as calculated from 
Eq. ([5} as a function of fc/fcp for Qcc = 0.5. The data labeled 
by the dashed (red) line correspond to A = 10"^ in Eq. (|16|l . 
while the data labeled by the solid (blue) line correspond to 
A = 1. Top panel: A = 10. Bottom panel: A = 100. In both 
panels the vertical (black) line indicates the point fc = fcp and 
the right axis scale refers to the A = 1 data. 



the leading term in the small x' expansion of Vm goes 
like x'™ . As we explain below, this m-dependence is re- 
sponsible for the main qualitative difference between the 
physics of quasiparticle mass renormalization in single- 
layer and bilayer graphene. 

For A = fcniax/^F large and A not too small, 
the quasiparticle mass renormalization is due entirely 
to the pseudospin dependent part of the self-energy. 
(The pseudospin-independent contribution, illustrated in 
Fig. [21 varies as k^ at small fc.) Neglecting Bq wc find 
that the renormalized velocity at fc = fcp is 



V* _ 1 dT.{k) 
V V dk 



k—k-c 
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= 1 



1 dfAjx) 

2"'='= dx 



(19) 



where 
/a (a:) 



dx'x'Vi{x, x') 



^ dx'x'^ 2F,{3/2, 1/2, 2, x'Vx^) 



rfx'- 2^1(3/2, 1/2, 2, xVa;'^) 



(20) 



It follows that 
dfA{x) 



dx 



2 7i 



dx'- 2^1(3/2, 1/2, 2, l/a;'2) 



1 



dx'^ 2^1(5/2,3/2, 3, l/x'2) . 
1 ^ 

(21) 



Because 



lim 2i^i(3/2,l/2,2, l/a;'2) 

x' — »-oo 

= lim 2-Fi(5/2,3/2,3,l/a;'2) = 1 , (22) 

x' — >oo 

we immediately see how the first term in Eq. (j2ip diverges 
logarithmically for large A: we find, to leading order in 
A for A ^ 00 



dfAix) 



lim 

A^cxD ox 



1 



ln(A) 



(23) 



and thus we recover the well-known^ asymptotic result 
for the velocity enhancement 



lim - = 1 + ^ ln(A) 



(24) 



III. 



VELOCITY RENORMALIZATION IN 
BILAYER GRAPHENE 



Bilayer graphene has four atoms per unit cell and 
its vr-electrons are therefore described by a model with 
four-component spinors. In order to focus on the simil- 
iarities and differences between single-layer and bilayer 
velocity renormalization we use a commonly adopted 
two-component model for a bilayer which applies^" at 
energies below the inter-layer tunneling scale and ex- 
plains the bilayer's unusualii quantum Hall effect. The 
two-component model for bilayer graphene has also 
been used to calculate the bilayer compressibilityi^ and 
the static nonintcracting density-density linear-response 
functioni^. 

The key differences between the two-band models of 
single layer and bilayer graphene are: i) the band disper- 
sion in the bilayer case is quadratic with an effective mass 



m = t/{2v'^) (where t is the inter-layer tunneling and v 
is the Fermi velocity in an isolated monolayer), ii) the 
chirality is J = 2 rather than J = 1 for bilayer graphene 
(see below), and iii) the intra-layer [V^^'' = Wfc] and inter- 
layer \y'^^ = v'^^ exp(— fed)] Coulomb interactions are 
different in the bilayer case. In our discussion of bilayer 
graphene we use a Thomas-Fermi-like potential for as 
in Eq. ()16p . to avoid the well-known mean-field-theory 
artifacts in Coulombic systems. We also employ a cut-off 
^max for the bilayer case, although we shall see that its 
role is less essential. (The two-band continuum model ul- 
traviolet cutoff scale in bilayer graphene is set by t which 
is smaller than the 7r-band width scale ^ v/a appropri- 
ate in the single-layer case.) Our calculations for bilayers 
have screening (A) and cut-off (/cmax/^F = A) parameters 
as in the single layer case, and are in addition depen- 
dent on the dimensionless inter-layer distance parameter 



dkr, 



which has a valued « 0.2. The Thomas- Fermi 



screening vector for bilayer graphene is given by gxF = 
me^/e so that gxr/^F = [i/(2f fcinax)]acoA = toeeA. Us- 
ing d = 3.35 A and t — 0.3 eV, we find that 

^ = toccA ~ 0.38 aocA . (25) 
fcp 

The mean-field theory calculations for the two-band 
model of bilayer graphene follow precisely the same lines 
as in the single-layer case. Eqs. ([S]) and ^ become: 

So(fe) = - / ^vi%U{k') (26) 



(27r)2 



and 



Zm 

d^k' 



(27) 



In Eq. (|27|l uj{k) = {cos{Jipk),sin{J(pk)) specifies the 
fc-dependence of the direction in pseudospin-space of 
the band-structure contribution to the effective mag- 
netic field. In the single-layer case the chirality J has 
the value J = 1 and Wi(fc) — k. It follows that only 
the pseudospin-dependent part of the Hartree-Fock self- 
energy differs between single-layer and bilayer cases: 



dk' 
2^ 



fc'/-(fc')^2^°^(fc,fc') 



(28) 



Because the two low-energy sites which comprise the bi- 
layer's pseudospin are in different layers, the pseudospin- 
dependent self-energy is proportional to the interlayer 
interaction V^^'^ . The change from the m = 1 angular 
component of the interaction in the single-layer case to 
the m = 2 component in the bilayer case follows imme- 
diately from the J = 1 to J = 2 band-chirality change. 
In Eq. 



vi'^HKk') = 



27r 



d9 



27r 
X exp(— 



{y/k^ + k'^ -2kk' cos{0)) 



2kk'co8{e)) . (29) 
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Using Eq. ((T3|) we can rewrite Eq. pS)) as 



(32) 



E(fc) 



1 

47r 



(30) 



In dimcnsionless units we have (scaling energies with vk-p, 
which is not the Fermi energy!) 



S(x) 



1 
2 



dx'a;'t/2*°^(:E,x') 



(31) 



where, as in the single-layer case, all wavevectors have 
been rescaled with fcp, i.e. x — k/k-p, x' = fe'/fep, and 
where we have introduced the dimensionless interaction 
vj^\ which is Kn^' measured in units of 2TTe^ / (ekp). 

The dependence of S(fc) on k/kp for the bilayer case is 
illustrated in Fig.[3]for acc = 0.5 and A = 10 and 100. In 
this plot we have shown S(fc) for two different values of 
the screening parameter A defined in Eq. p^ : A = 10^^, 
corresponding to an essentially unscreened Coulomb po- 
tential, and A = 1 corresponding to the Thomas-Fermi 
screened potential. The small k behavior is not strongly 
influenced by screening and can be understood analyti- 
cally. For large A we can expand the exponential that 
enters in Eq. ((29l) in powers of d/A: 



exp 



A 



^yx'^ + - 2xx' cos{9)^ 



^ 1 — —y/x'^+x''^ — 2xx' cos{6) 

+ i (j^ [x^ + x'^ - 2xx' cos{0)] + 0{{d/Af) 



For unscreened Coulomb interactions 



V(^\x,x') 



dO 
— ( 

2tt 



X exp 



\J x^ + x'"^ — 2xx' cos{6) 
^- ^Jx'^ + a;'2 - 2xx' cos(0) 



(33) 



This implies that the expansion of vS^\x,x') starts at 
second order in c?/A for m > 1: the term proportional 
to -y/x^ + — 2xx' cos(6') , which is first order in d/A, is 
canceled by the denominator in the integrand in Eq. (j33|) . 
Thus the angular average in Eq. p3p of this term would 
give us a finite result only for m — Q. 

Disgarding this term, we obtain the following expan- 
sion of V^k'^x^x') in powers of d/A valid for to > 1 : 



1 + 5 (|) (- + ^") 



Vm{x,x') 



^XX' [V„i+i{x, x') + V,n^i{x, x')] 



(34) 



Here the coefficients Vmix, x') are given by Eq. pB]). Us- 
ing this expression in Eq. (|3ip we find 



E(x) 
vkp 



dx'x'V2{x, x') + -etc 



A 



2 „A 



x^ I dx'x'V2ix,x') 

2 „A 



A 



dx'x'^V2{x,x')- -a^J-] X dx'x'^[Vi{x,x') + V3{x,x')] 



(35) 



To understand the limit a; ^ of the bilayer Hartree- 
Fock self-energy in Eq. ([35|) we use that 



V,n{x^O,x') = 



X 



1 • 3 • 5 • . . . • (2m- 1) 



f,/7n-\-l 



2™to! 



(36) 



Inserting this expression in Eq. (I35|) . carrying out the 
integrations over x' , and neglecting terms that go to zero 
faster than x'^ for x — > 0, we finally find 



S](a; ^ 0) 
vkp 



A- 1 
A 



16 

l-a,^x'd' ^+0{x^) . (37) 



doping limit), while the second term which contains the 
layer separation dependence goes to zero. As in the sin- 
gle layer case, the dominant contribution to the exchange 
self-energy at small wavevectors originates from interac- 
tions with states deep in the negative energy sea which 
are not sensitive to screening. 

The renormalized mass to* in bilayer graphene is de- 
fined via 



kp 

TO* 



d_ 

dk 

kp 

TO 



2to 
1 
2 



Bo{k) + S(fc) 
dTAix) 



k—kp 



dx 



(38) 



Note that the first term is finite in the limit A — > cx) (zero Restoring the pseudospin-independent self-energy contri- 
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FIG. 4: (Color online) The Hartree-Fock renormalized mass 
m* /m for bilayer graphene [calculated from Eq. (|40|l ] as a 
function of Qco for A — 10~^ and various values of A. From top 
to bottom, A = 10 (squares), 20 (circles), 50 (triangles), and 
100 (hexagons). Note how interactions suppress the quasipar- 
ticle effective mass. 



enters the Thomas-Fermi screening vector. We plot the 
ratio m* /m obtained from this expression as a function 
of c^ee for A = 10~^ in Fig. [4] while in Fig. [5] we plot 
m* /m as a function of density for various values of the 
screening parameter A. In this last plot we have taken 
into account the spin and valley degeneracy factor g = 4 
of graphene by letting gxF giTF- 



FIG. 3: (Color online) The pseudospin-dependent part of the 
Hartree-Fock self-energy S(fc) (in units of vkv) for bilayer 
graphene as a function of fc/fep (for wavevectors up to the cut- 
off) for dec = 0.5. The data labeled by the dashed (red) line 
correspond to A = 10~^ in Eq. (|16|l . while the data labeled by 
the solid (blue) line correspond to A = 1. Top panel: A = 10. 
The vertical (black) line indicates the point k = kw- Bottom 
panel: A = 100. 



bution, which is identical in single layer and bilayer cases, 
we obtain that 

TAix) = / dx'x'V^^\x,x') -2 dx' x'Vo{x,x') 



dx' x'Vo{x, x') . 



It follows that 



m 
m 



1 



1 dV^jx) 



(39) 



(40) 



where t, defined above Eq. (l25l) . is the same quantity that 



IV. SUMMARY AND DISCUSSION 

In summary, we have presented a comparative study 
of Fermi surface effective mass renormalization in single- 
layer and AB-stacked bilayer graphene. Interband ex- 
change adds a new twist to electron-electron-interaction 
velocity renormalization physics compared to ordinary 
two-dimensional electron systems^ii^, one which always 
leads to strongly enhanced velocities and correspond- 
ingly suppressed masses. We have shown analytically [see 
Eq. (|37p ] and numerically (see Fig. [5l especially the in- 
set) that although this mechanism causes the mass to 
approach zero at low carrier densities (A ^ cxd) in single- 
layer graphene, the renormalized bilayer mass remains 
finite. The two cases are different because of the differ- 
ence in the pseudospin-chirality of their quasiparticles, 
J = 1 in the single layer case and J = 2 in the bilayer 
case. Although the bilayer mass is finite, it is still esti- 
mated to be very strongly suppressed, depending on the 
approximation used for screening. 

This behavior should be contrasted with what happens 
in the case of conventional parabolic-band electrons liq- 
uids in Si-MOSFETsi^'l^ or in GaAsil, where the quasi- 
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2 4 6 8 10 12 
n (in units of 10-'-^ cm^^) 

FIG. 5: (Color online) The Hartree-Fock renormalized mass 
m* /m for bilayer graphene [calculated from Eq. (|40p ] as a 
function of density n (in units of lO^'^ cm~^) for dee = 0.5. 
The (blue) circles correspond to A = 1, the (red) squares to 
A = 10~^, while the (green) triangles to A = 10"*. In the inset 
we have reported logjQ(m*/m) as a function of logj^Q(n) for 
A = 10~^ and A = 10~*. Note that the bilayer renormalized 
mass always saturates at a finite value for n ^ 0. 



particle velocity is generally suppressed by interactions 
(except in a narrow interval of very high densities - small 
values of the Wigner-Seitz density parameter - which 
depends in a complicated way on extrinsic effects such as 
the quantum- well finite thicknessi^) . The true value for 
m* as a function of density for graphene bilayers proba- 
bly occurs for a value of the screening parameter A ~ 0.5, 
which represents a compromise between using Thomas- 
Fermi theory and complete neglecting screening effects. 
It should be emphasized however that a reliable numer- 
ical estimate for m* as a function of density in bilayers 



lies outside the scope of this work and, as in the ordinary 
two-dimensional electron gasF^i^^^, may ultimately re- 
quire elaborate calculations that are informed by quan- 
tum Monte Carlo numerical studies. 

We conclude by commenting briefly on the possibil- 
ity of broken symmetry states in single-layer and bi- 
layer graphene systems. Enhanced quasiparticle veloc- 
ities tend to increase the kinetic energy cost of mov- 
ing the chemical potential further away from the Dirac 
point. The same physics which leads to velocity en- 
hancement therefore also tends to suppress spontaneous 
spin-polarization, the type of instability most often con- 
templated^ii^ in ordinary low-carrier-density electron gas 
systems. Because of the pseudospin structure of single- 
layer and bilayer graphene, interactions instead favor 
chiral-symmetry-breaking-type scenarios, as we have ex- 
plained elsewhere^. Chiral symmetry breaking, which 
would lead to a spontaneous gap at the Dirac point and 
in the bilayer case to spontaneous charge transfer be- 
tween layers^, is most likely to occur in bilayer graphene. 
If these instabilities occur, they are confined to a region 
of low carrier density-^ and the discussion of the current 
paper would apply only outside of this regime. In bi- 
layer graphene the occurrence or absence of chiral sym- 
metry breaking may depend on higher neighbor inter- 
layer hopping processes not included in the model dis- 
cussed here. Recent lattice Monte Carlo calculations by 
Drut and Lahde^l suggest that these type of instabilities 
even occur for the J 1 single-layer case. 
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